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Abstract 

Cosmic Microwave Background (CMB) anisotropies that result from quantum fluctuations during 
inflation are explored and the impact of their "cosmic variance" on the ability to use existing data to 
probe inflationary models is studied. We calculate the rms temperature fluctuation, and its cosmic 
variance, for a number of experiments and for models with primordial power spectra which range 
from n = | to 1. We find: (1) cosmic variance obscures the information which can be extracted, 
so a comparison of the rms temperature fluctuation on small scales with the COBE result can fix 
n to only « ±0.2 at best; (2) measurements of the rms fluctuation on 1° scales may not allow one 
to unambiguously infer the tensor contribution to the COBE anisotropy; (3) comparison of this 
contribution with the predictions of inflation are ambiguous if the quadrupole anisotropy alone is 
utilized. We discuss means for minimizing the uncertainty due to cosmic variance in comparisons 
between experiments. 
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1. Introduction 

One of the central purposes of measuring CMB anisotropies is to gain information on the power 
spectrum of primordial density fluctuations, and from this, information on the physical processes 
in the early universe which may have produced such a spectrum. The well studied mechanism 
of inflation allows an a priori calculation of the spectrum of fluctuations. Both scalar and tensor 
fluctuations arise from calculable quantum fluctuations in elementary fields during an inflationary 
phase in the early universe, and any previously existing density fluctuations are inflated away. We 
explore here what constraints on the primordial power spectrum can be obtained, both in princi- 
ple and in practice, from the COBE DMR observations in combination with observations of CMB 
anisotropies on scales as small as 1°. Comparison with large scale structure analyses (e.g. 5T/T\sw 
as inferred from POTENT (Bertschinger & Dekel 1989; Bertschinger, Gorski & Dekel 1990) ), could 
provide useful additional constraints although at present there are large theoretical and observa- 
tional uncertainties. 

When comparing theoretical predictions with observations, particular attention should be paid 
to the fact that CMB anisotropies resulting from inflation are stochastic. Because of the quantum 
nature of the process by which fluctuations are generated, only observations over an ensemble of 
universes can allow the unambiguous measurement of fundamental microphysical parameters, even 
in principle^]. Constrained to observe only one universe, there remains an irremovable uncertainty in 
our ability to relate certain CMB measurements, no matter how precise, to predictions of inflationary 
models (such as the normalization and slope of the power spectrum). For example, inflation predicts 
a probability distribution for ratios of the moments of CMB temperature fluctuations (or functions 
thereof), and not the ratios themselves. This has important consequences for the normalization of 
the power spectrum. We term the induced uncertainty the "cosmic variance", not to be confused 
with other usages of this term. Among our results will be a description of the effect of this variance 
on the constraints which can be derived from present and future observations of the CMB. 

2. The Primordial Power Spectrum and Inflation 

The remarkable observations of primordial anisotropies in the CMB by the DMR instrument 
aboard COBE (Smoot et al. 1992) offered cosmologists the first "unprocessed" glimpse of what 
may be the initial conditions for the present observed Friedman-Robertson- Walker expansion of the 
universe. Inferring the primordial power spectrum of fluctuations based on observations of large 
scale structure had previously required detailed assumptions about both the equation of state of 
matter inside the horizon and the relationship between the observed structure and the underlying 
distribution of mass in the universe. Now model predictions can be normalized at COBE scales, 
which represent modes which had not crossed the horizon at recombination and hence could not have 
been processed by causal microphysical factors. This represents a significant theoretical advance. 

2 We shall discuss this issue in more detail later in the text. 
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Moreover, the scales observed by COBE in combination with small scale observations provide a 
large lever arm for exploring the scale dependence of the power spectrum. 

A central quantity of interest has been the slope of the power spectrum of fluctuations. Defin- 
ing the power spectrum, P(k) in terms of the Fourier transform of the density fluctuation: P(k) = 
\J(Sp/ p) exp(ik ■ x)d 3 x\ 2 , one finds that the rms density fluctuations on scales with comoving wavenum- 
ber, k can be written (e.g. Peebles 1980) 



One conventionally considers P(k) oc k n . This assumes a scale-free spectrum, namely that there 
is no preferred scale for fluctuations set by the underlying physics. Further, n = 1 represents a scale- 
invariant spectrum as predicted in standard inflation (Steinhardt & Turner 1984), up to logarithmic 
corrections. It has been stressed, coincident with COBE, that a number of inflationary models 
which attempt to circumvent the problems of new inflation also predict measurable deviations from 
a scale invariant spectrum, i.e they can predict n < 1 (Davis et al. 1992; Gelb et al. 1992; Liddle & 
Lyth 1992a,b; Lidsey & Coles 1992; Lucchin, Matarrese & Mollerach 1992; Salopek 1992). COBE 
reports a limit on this spectral index of n = 1.1 ± 0.5, which is consistent with scale invariance as 
well as with significant deviations from scale invariance. Clearly to get an optimum handle on this 
spectrum one would like to compare large-scale and small-scale anisotropies in the CMB. 

3. CMB anisotropies 

It is useful to decompose the temperature fluctuations in spherical harmonics AT/T(6, <p) = 
Y J im a imYim{9,4 > ) an d to define the rotationally symmetric quantity of = J2 m \ a im\ 2 - We choose to 
consider here the so-called rms fluctuation, or correlation function at "zero-lag" , as our contact with 
experiment. This has a simple interpretation and is sensitive to a range of physics from both large 
and small scales. To be precise we consider 



where the dipole has been subtracted and W\ is a "window function", which incorporates the 
experimental sensitivity to each moment. The window functions appropriate to the experiments 
that we will consider in this paper are described, and derived, in appendix A. As written, AT/T 
has contributions (aj)| from scalar (density fluctuations) and {ai)\ from tensor (gravity wave) 
modes. For the COBE experiment, the window function is simply due to the gaussian smoothing, 
Wi = exp [—(I + 1/2)V 2 ] with a = 0.425 x 10° the smoothing scale (the factor 0.425 relates the 
gaussian width a to the FWHM quoted by the COBE group). We will first concentrate on this 
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"window", where AT/T then represents the rms temperature fluctuation over the observed sky. 
Understanding the cosmic variance of the COBE measurement is important if COBE is to be used 
to normalize the power spectrum to compare with other experiments. 

Inflation predicts that the a\ m are independent gaussian random variables (=>- af is x 2 -distributed 
with 21 + 1 d.o.f.) so the rms fluctuation as defined is itself a stochastic quantity (Abbott & 
Wise 1984; Krauss & White 1992; White 1992). It might be thought that measurements made over 
the entire sky afford the opportunity to make independent measurements of the quantities which 
have a cosmic variance, and so reduce it. However from the above, it is clear that this cannot be the 
case. No measurement of AT/T on our sky, no matter how accurate nor how many independent 
measurements it corresponds to, can reduce the cosmic variance for a given /. Determining each af 
requires an all-sky average (we will not address the important issue of the effects of incomplete sky 
coverage in this paper as we are interested here in the irremovable cosmic variance. The issue of 
finite sky coverage is addressed in Scott et al. 1993.). The error in measuring af will be reduced by 
more sky measurements, but not the uncertainty in going from af to the average over an ensemble of 
universes {af) ens . (However, as we will see, window functions relevant to small scale measurements 
can give dominant contributions to AT/T from higher modes, which have more degrees of freedom 
and smaller cosmic variance, so that the effect of this uncertainty can in principle be lessened, 
provided one has adequate sky coverage.) 

It is not possible in general to find a closed form analytic expression for the distribution of 
AT/T. We outline a simple and accurate numerical procedure for obtaining the distribution in 
appendix B. It is straightforward, however, to calculate the first few moments of the distribution 
(since AT/T is a sum of independent, x 2 -distributed terms) which give us an estimate of the size 
of the variance. Write 

/AT\ 2 of _ 2 — (of) 

\ T J Y ( a i f) 2/ + 1 4tt v ; 

so (AT/T) 2 = ((AT/T) 2 ) = £j(Z+l/2)A*. If we define AW = (((AT/T) 2 -(AT/T) 2 ) n ) then we can 
use A( 2 ) as an estimate of the (cosmic) variance of AT/T and A® as an indicator of the asymmetry 
of the distribution about the mean (and departures from normality). It is straightforward to 
compute^ 

A (2) = + V2)A 2 A (3) - 2J2(1 + l/2)Af . (4) 

i i 

For n = 1 Sachs- Wolfe fluctuations ((af) ~ (/+1/2)" 1 , 1 > 1), the rms fluctuation (AT /T(a)) 2 ~ 

3 Our result for A^ 2 ) differs from Bond & Efstathiou (1987) by a factor of 2. 
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log a while the higher moments are all finite^. Thus the "cosmic variance" becomes less important 
at smaller scales. It is at these small angular scales (less than about 2°) however that the theoreti- 
cal uncertainties due to cosmological model dependence and non-linear effects become progressively 
more severe. 

To better illustrate the form of the distribution for AT/T we have compared two approximations 
to the exact distribution (which we checked by doing a Monte-Carlo analysis in which the af 
were drawn at random from X21+1 distributions and (AT/T) 2 computed using equation [|). In 
the case of only a scalar component in fluctuations of the n = 1 Sachs- Wolfe form (specifically 
(af)s = (/ + 1/2) -1 for I > 2), we show in Figure la the distribution of AT/T smoothed to 10° 
obtained from the method outlined in appendix B and the Monte-Carlo calculation (see also Table 
H). These are in good agreement. Also shown are fits obtained from matching the mean and 
variance of the distribution to a x 2 form {P{x) oc x Cl e~ C2X ) and a gaussian. As can be seen, it is 
not well fit by a gaussian. The fit to a single x 2 form gets progressively better for larger angles, or 
tilted spectra. This is because the sum is dominated by the lower multipoles in these cases and is 
therefore closest to being represented by a single x 2 distribution. It is important to note that if one 
attempts to fix the normalization of the power spectra by fitting AT/T to the COBE measurement, 
there is an uncertainty due to cosmic variance of ~ 20% (figure la). This normalization uncertainty 
will affect any comparison with smaller scale experiments and reduce what can be inferred. In the 
future it may be better to normalize the spectrum not by the COBE AT/T but by something less 
sensitive to the low I multipoles. We will return to this point later. 

For smoothing angles smaller than 10°, one can use the second moment of the distribution as 
an equivalent gaussian error for the purposes of computing confidence levels. We again stress that 
the distribution is not gaussian but for quantities depending only on the asymptotic form of the 
distribution (e.g. 90, 95, 99% confidence regions) such an approximation is adequate. Nevertheless 
this skewness should not be ignored when utilizing measurements of either AT/T or the correlation 
function C(6) to place limits on microphysical parameters such as the scale of inflation. 

4. The Anisotropy at 1° Scales 

On degree scales, for the scalar fluctuations, one must consider in addition to the Sachs- Wolfe 
terms, other terms such as the Doppler contribution (from motions on the last scattering surface). 
These contributions depend on the cosmological model (e.g. the Hubble constant, VLb and the re- 
ionization history of the universe), which governs the growth of fluctuations inside the horizon, and 
generally require numerical solution of the Boltzmann equation for the evolving photon distribution 
for an accurate estimation. To calculate (of) we have used the radiation "power spectrum" k 3 W^(k) 
of Vittorio & Silk (1992), which includes a "transfer function" for the evolution of modes which have 

In practice of course the rms fluctuation will not diverge as a — > 0. Physics associated with the size and width 
of the last scattering surface damps the af for high I. This can also be aided by spectral tilt. 
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crossed the horizon. As an illustrative example of the non-Sachs- Wolfe terms, we calculated the 
(af) for a CDM model with Q = 1, h = 0.5 and Qb = 0.03, 0.1 (consistent with the range favoured 
by nucleosynthesis: Smith et al. 1992, Krauss & Romanelli 1990, Walker et al. 1991). For lwe 
multiplied the Sachs- Wolfe result by the ratio of the full calculation to Sachs- Wolfe for n — 1. This 
approximation is expected to be good to better than 10%, as can be seen from analytic work (Atrio- 
Barandela & Doroshkevich 1993). We normalized the af to fit the COBE central value at 10°. The 
power spectra are shown in Figure 2 for n = 1.0 for the range of k = IHq/2 of relevance here. (An 
unprocessed Sachs- Wolfe power spectrum would go as k 71 " 1 .) The (af) are an integral of the "power 
spectrum" with measure jf(2k/H Q )dk/k — in effect an average of the "power spectrum" around 
k pa IHq/2. Hence the 'Doppler enhancement' of the power spectrum on degree scales translates to 
an increase in the (af) for large I over the Sachs- Wolfe result. To obtain ((AT/T) 2 ), one performs 
a weighted sum of (af) with weight W\ as in equation (f2|). Plotted in Figure 2 are the window 
functions for the COBE, MIT, UCSB South Pole (SP91) and MAX experiments to show the range 
of multipoles probed by eacbQ. 

In Figure 3 we display the predicted value of AT/T (for Q B = 0.03,0.10) for SP91 vs n, for a 
spectrum which is normalized to fit the COBE measured AT/T. (For a more detailed comparison of 
SP91 with COBE in the context of standard CDM see Gorski, Stompor & Juszkiewicz 1992). The 
solid lines represent the 90% confidence level predictions assuming a tensor component is present 
in the amount predicted by power-law inflation (Davis et al. 1993) and the dashed lines assume 
that the tensor component is absent (see section 5). Note that the entire uncertainty in the values 
shown comes from cosmic variance of the COBE result, with present observational uncertainty 
ignored and no uncertainty for SP91. This curve allows a determination of the constraints on n 
which may be possible from a comparison of SP91 and COBE. Concentrating for now on the solid 
curves, giving the full result (including tensor contributions to be described next), we see that even 
an exact measurement of AT/T by SP91 may only constrain n to within pa ±0.2, due to the COBE 
cosmic variance. Including a SP91 measurement uncertainty might, depending on its value and the 
measured AT/T, allow the whole range n = 0.5 — 1. The corresponding predictions for the MAX 
experiment are shown in Figure 4. 

In principle one might reduce the uncertainty on the normalization of the (af) inferred from 
COBE (and hence the theoretical uncertainty on smaller scales) by using not just AT/T but the 
whole correlation function C(9), which contains more information. Because C(9) oc J2i afPi{cos 6)Wi 
probes a different linear combination of af at each 9, one might hope for increased sensitivity to 
higher I modes which have smaller cosmic variance. Nevertheless, in practice the present sensitivity 
is such that most of the signal comes from small /, so that any treatment of the existing data cannot 
decrease our estimate of the errors to the point where cosmic variance is irrelevant. The relative 

5 Recall in our convention the window function is convolved with af which fall with I in the Sachs- Wolfe region. 
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uncertainty and size of the (af) are shown in Table [I] for reference. 

As Figure 2 shows, the MIT experiment probes a larger range of I than COBE and the effects 
of cosmic variance on the normalization inferred from this experiment are less. We find that the 
normalization uncertainty is ~ 13% (thus the MIT experiment would be more useful than COBE 
for fixing the normalization of the spectrum if this experiment were to have extensive sky coverage 
and reduced experimental errors). As shown in Figure lb (AT/T) 2 still has significant skewness 
and is not well approximated by a gaussian of mean ((AT/T) 2 ). Again the correlation function 
(assuming small errors at each angle) would contain more information than AT/T, but would be 
more difficult to analyze. Even so it would be useful to have a measurement of the anisotropy on 
scales large enough that sensitivities to the doppler peak (Qb and amount of tensor component) 
and re-ionization are small, while being such that low / modes to not contribute significantly e.g. 
a "full sky MIT" correlation function or an experiment intermediate between Tenerife (Davies et 
al. 1987), which measures fluctuations on 8° scales, and SP91 on degree scales. 

5. Tensor versus Scalar modes 

It was pointed out in Krauss & White (1992) that gravitational wave (i.e tensor modes) could 
play an important role in interpreting the COBE anisotropy if it stems from inflation, so that 
these should not be ignored. As a result, several groups examined the ratio of scalar to tensor 
modes predicted to result from various inflationary models (Davis et al. 1992; Dolgov & Silk 1993; 
Krauss 1992; Salopek 1992). In Davis et al. (1992), a general connection was pointed out between 
the scalar to tensor ratio predicted by inflation and the spectral tilt as given by n ^ 1. For power 
law inflation this relation takes the approximate form n pa 1 — 1/7(T/S) (Davis et al. 1992) where 
{T/S) represents the ratio of the expectation values of the tensor and scalar quadrupole moments. 
As a result it appears that this ratio, if it were measurable, could provide important information 
about an inflationary epoch in the early universe (Davis et al. 1992). It is important to note however 
that while the ratio of the expectation values are fixed by the equation of state during inflation, 
the separate cosmic variance can alter the actual "observed" ratio of quadrupoles significantly from 
that predicted by T/S. For two x 2 -distributed random variables with 21 + 1 d.o.f., the ratio has the 
distribution 

P(R)dR- F(2/ + 1) Rl ~ 1/2 B -A1±_ (,) 
[ ' ~ T 2 (/ + 1/2) (1 + i?) 2/+1 K ~ (a 2 )/(af) 1 } 

In Figure 5 we have plotted P(R) and the integral J R P(r)dr for the quadrupole (1 = 2). Note that 
the distribution is peaked at a value below 1 but has a significant large R tail. With P(R\ < R < 
R2) = 95% we find that R can range from 0.2 — 5! Hence fluctuations with [a^r/ '{ a V)s 5 times 
(or 1/5) the ratio (a^r/ '(cfys are not improbable. Thus we see that a direct measurement of the 
quadrupole ratio (a|)T / '(al) s could not provide much information on the actual equation of state 
during inflation! One needs to look at quantities (like AT/T) that have smaller cosmic variance 
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such as the power spectrum normalized quadrupoles Qrms-ps (Smoot et al. 1992). In this case the 
conversion to T/S is much less uncertain. The term quadrupole in this context should therefore be 
interpreted with care. 

Nevertheless any confirmation of a tensor component to the CMB anisotropy on COBE scales or 
below would be of great significance for the interpretation of the COBE measurement, whether or 
not it pins down models. However, our results point out the difficulty in making such a confirmation 
based on comparing CMB anisotropics on 1° scales with the COBE result. The solid curve in Figure 
3 includes the tensor contribution in the amount prescribed by power law inflation as a function of 
n. If the tensor contribution were absent for some reason (dashed curve), the predicted AT/T at 
small angles would be increased. This is because the (af)r begin to damp (to be conservative) at 
I ~ \/l + Zd e c ~ 30 due to redshifting of the gravity waves (Starobinsky 1985; Turner et al. 1993). 
It is therefore smaller than the scalar for high I modes, so that its relative contribution to COBE 
is larger than it is to small angle experiments [Note that the spectral tilt required to have a large 
tensor component increases the contribution of low / modes to the sum, reducing the sensitivity 
to the Doppler peak and makes tensor and scalar contributions similar.] While this provides some 
sensitivity to a tensor component in principle, for SP91 the maximum variation is only slightly 
greater than the (90% confidence level) irremovable uncertainty due to cosmic variance. Including 
the uncertainty in VLb (and n) it seems that a comparison of 1° scale measurements with COBE can 
at best marginally probe the expected tensor contribution to AT/T. The situation on still smaller 
angular scales is less pessimistic in this regard, but one should remember that these scales are most 
sensitive to non-linear effects and assumptions about the cosmological model, as shown in Figure 4. 

The results of this section, and those of section 4, demonstrate the essential problem involved 
in comparing the COBE results to those of small angular scale experiments. If the window func- 
tion of such experiments probes different / than COBE, the 20% COBE cosmic variance can blur 
the comparison. To overcome this normalization uncertainty, we suggest using the COBE correla- 
tion function or AT/T as measured by the MIT experiment to fix (a|). In the latter the window 
functions for small scale experiments can overlap the MIT window function significantly, and corre- 
lations between the variance affecting the normalization and the variance in the small angular scale 
predictions must be taken into account. 

6. Cosmological Model Dependence 

As alluded to earlier, while the COBE measurement probes scales outside the horizon today 
and is mostly insensitive to details of the cosmological model, on smaller scales the predictions 
are model-dependent. Apart from uncertainties in cosmological parameters (such as Hq and Qb) 
which induce uncertainties in the prediction of AT/T, the major unknown arises from the possiblity 
of early re-ionization and re-scattering of the primoridial CMB fluctuations. The unbiased CDM 
model, precisely that normalized to COBE, allows adequate non-linearity at early epochs that re- 
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ionization and generation of a secondary "last scattering surface" (LSS) is not implausible (Tegmark 
& Silk 1993). The ionization of hydrogen imposes such modest energetic requirements that rare 
fluctuations of low mass can provide a plausible source of ionizing photons or particles (an effect that 
must have occurred by z = 5 to satisfy the Gunn- Peterson constraint). Such re-ionization would 
move the LSS in to lower redshift and reprocess fluctuations on scales smaller than the horizon 
size at the new LSS. Indeed for a range of (high) CDM normalizations with n = 1 one can reduce 
degree scale fluctuations by a factor ~ 2 — 3. However if the fluctuations are small (corresponding 
to biassing factors b ~ 2), the spectrum is tilted or is that of a mixed dark matter model, it appears 
that degree scale fluctuations remain largely unaffected. 

Secondary fluctuations induced on very small (arc-minute) scales and their detection could even- 
tually provide a useful constraint on the last scattering surface. Until then however primary fluc- 
tuations on degree scales are potentially more sensitive to the last scattering surface (re-ionization) 
than to the possible presence of tensor modes. 

7. Conclusions 

After the initial flurry of excitement over the COBE data and its potential ability — in com- 
bination with other CMB measurements — to probe the details of models for the generation of 
primordial density fluctuations, it is time to realistically assess the practical and theoretical limita- 
tions inherent in any such procedure. This is particularly important as new measurements of CMB 
anisotropies on 1° scales are expected. Our results indicate that the statistics of CMB anisotropies 
predicted by inflation play a fundamental and irremovable role in limiting the information we can 
gain from the data. Cosmic variance, even with an exact COBE result, still provides uncertainty in 
the extracted spectral index n. By incorporating this variance, we have been able to quantitatively 
display how limits on AT/T at small angles can translate to limits on n. Going to smaller scales, 
such as those comparable to that of observed large scale structures, would provide a longer lever 
arm and in principle tighter constraints on n. However the theoretical uncertainties due to causal 
effects and their relation to the growth of fluctuations are much larger in this case. Finally, while 
observation of a tensor component to the CMB anisotropy would be of great interest, we have shown 
that the predicted levels are not easily extractible from the data. 

To reduce the effects of cosmic variance we suggest that the full correlation function of COBE or 
the MIT experiment would give a theoretically cleaner normalization of the power spectrum than 
the COBE AT/T. This can then be used with the South Pole or MAX experiments to constrain 
models of structure formation or cosmology. The normalization uncertainty in this case can be 
reduced by perhaps a factor of 2. However this would still make a detection of a gravitational wave 
signal very difficult, coupled as it is with the uncertainties in n, Qb and re-ionization history. 

Inflation is an attractive theory in part because it makes unambiguous predictions which can 
be tested. Similarly, measurements of the CMB anisotropy still provide the best hope we have of 
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extracting information on the primordial power spectrum of density perturbations. As we have 
shown however, definitive answers will require careful attention to the theoretical uncertainties. 
Inflation predicts probability distributions for observed CMB anisotropies. It will be important to 
determine what CMB measurements can rule out conclusively. On the other hand, it is equally 
important not to overstate the case. 
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Appendix A 

The sensitivity of experiments to the spherical harmonic decomposition of the CMB tempera- 
ture fluctuations is conventionally described in terms of window functions Wi as in equation (g). 
For completeness, and to fix our conventions, we describe in this appendix the window functions 
appropriate to the COBE, MIT, UCSB South Pole and MAX experiments. The functions are shown 
in Figure 2. 

COBE and MIT 

A simple way to view the window functions is in terms of a mapping Y\ m — > Ti m Yi m (no sum) 
describing the weighting the experiment gives to each mode. An example makes this clear. The 
most straightforward such mapping is that caused by a gaussian smoothing function of width a <C 1 
(which could be due to antenna response or applied to the data "by hand") 



Yim{0o, 0o) -> / 2 2 exp 



-0 2 



2a 2 



Y lm (tt + n) 



exp 



-(Z + 1/2)V 



Y 



lm{V0, <P0) 



2 

(for details see White, 1992). If we substitute 5T/T = J2i m a im T imYim into 

where (■ ■ •) indicates an average over the whole sky, and use the orthonormality of the Y\ m we 
obtain the window for the COBE experiment (Smoot et al. 1992) W\ = exp[— (I + l/2) 2 a 2 ] shown in 
Figure 2. In this case AT/T corresponds to an rms temperature fluctuation over the observed sky, 
smoothed on a scale a = 0.425 x 10°. (The COBE maps with no additional smoothing correspond 
to a — 0.425 x 7°, not shown in Figure 2). We find the same factor appearing in the other window 
functions due to the finite resolution of the antenna (which we assume to be gaussian). In what 
follows we will make the approximation of a "perfect" antenna, knowing we can put the antenna 
response in at the end. 

The MIT experiment (Page et al. 1990) also calculates AT/T from a map of the sky and so the 
form of the window is identical to that of COBE. The smoothing scale is a = 0.425 x 3.8°. 

South Pole and MAX 

Now consider the UCSB South Pole experiment (Gaier et al. 1991) and focus on the measurement 
of the so-called (AT/T) rms on one patch of the sky, centered at 9q, (fio- (Since it is our purpose here 
to show the effect of cosmic variance on inferences made from this experiment, and not to fit the 
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data, we will not be concerned with multiple, correlated patches of the sky. For a discussion of 
the fitting of the data see Bond et al. (1991); Dodelson & Jubas (1992). Note that in Bond et 
al. (1991) a factor of An/ (21 + 1) is missing from equation (1).) The temperature difference at 9 , O 
is measured by "chopping" the beam at fixed = 9 through <p(t) = <p + a sin(27rz/t) and sampling 
5T/T(6,<f)) with weight sign(0(t) — O )- Recalling the temperature difference is twice the weighted 
average, the mapping we want is then (see Dodelson & Jubas 1992) 

Yi m (0oAo) -> 2i/ [ 1/U dt sign(0(t) - fo)Y lm (e,<i>(t)) 



o 

= 2iH (ma)Yi m (9 , (f> ) 

where in the last line we have written the integral dt as an integral d((f>(t) — (fio), extracted the 
exp[im(f)} dependence from Yi m and used the identity (Gradshteyn & Ryzhik 1980) 

for the Struve function. At this point we could proceed to define the sky rms fluctuation by 
averaging over Q , O as before (now with an m dependent Tj m ). However the "rms" AT/T quoted 
by the SP group does not come from this procedure. Instead we now argue that the small angular 
scales probed by this experiment make it sensitive only to high / multipoles. These multipoles have 
a small cosmic variance, thus our universe at 9q, 4>q should be similar to the ensemble average of all 
universes at the same point. We can then define (AT/T) 2 = ((ST/T) 2 ) ens and use (dropping the 
subscript ens) 

( a l'm'0>lm) = (2Z +1) h' l ^ m ' m 

to write (AT/T) 2 = J2im( a f) VimYi m \ 2 / (21 + 1). Reinstating the finite beam width the window so 
defined is then 



' 47T 



21 ma] 



Wi= E ^ m = ^r— rexp[-(/ + l/2)V] £ |y /m (# o ,0o)| 2 4tf o 2 

m=— I V ' 1 m=—l 

For the configuration of (Gaier et al. 1991) O = 27.75°, asin# = 1-5° and a = 0.425 x 1.65°. 

If the effects of sky rotation and smooth vs. step scanning are ignored the form of the window 
function for the MAX experiment is very similar to SP91. However the MAX experiment uses a 
phase lock-in which extracts the frequency component of the sky signal which matches the chop 
frequency v. Extracting only this component corresponds to replacing the weight sign(0(t) — <f>o) 
with (<f)(t) — <fro)/a in the Y\ m mapping defined above. The remaining integral over dip is then not 
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Hq but a Bessel function: ±(7r/2) Ji(|ma|). The window function is as above for SP91 but with 
2Hq — > 7rJx- [This explains why our window function and predictions for MAX differ from that of 
other authors.] The parameters (for the /z-Pegasii scan, Meinhold et al. 1992) are 9 = 90°, a = 
0.65°, a = 0.425 x 0.5°. 

We can test our assumption that the cosmic variance is "small". Recalling that the ai m are 
gaussian independent complex variables it is straightforward to show that for an experiment covering 
the whole sky0 

' -AT_x 2 \_^(af). 




T 




I'm 



47T 



W, 



I'm 



where the Yi m factor is to be removed from W[ m . Notice that these reduce to equation when 
Wi m is m-independent. Keeping in mind that the distribution of AT/T will not in general be 
gaussian, this allows us to derive an estimate for the size of the cosmic variance associated with 
each experiment. We find the uncertainty (at IV") for an n = 1 spectrum to be 4 and 1% for the 
SP91 and MAX windows respectively. This would be larger for n < 1. 



Appendix B 

The probability distribution for the correlation function predicted by cosmic variance is not well 
fit by a gaussian or \ 2 form. In this appendix we derive a simple expression for the distribution 
which is amenable to numerical computation. We hope that this will facilitate comparisons of 
theory with data when they become available. 

For now let us consider the correlation function predicted by a one component spectrum (either 
tensor or scalar). We will keep the 9 dependence here for generality, and defer consideration of two 
component spectra till later. For COBE and MIT, by definition 

1 oo oo 

C{6) = — X>fP,(cos 9)W l = j>y, 

47F 1=2 1=2 

where we have defined yi = af/af with of = (af)/(l + 1/2) and q = (of /47r)Pz(cos 9)Wi. (Note q 
reduces to of equation @ for 9 = 0) For the Sachs- Wolfe part of the spectrum the q ~ l~ 2 W[ 

6 At the other extreme if the experiment looks at only one point on the sky the variance is twice the square of 
the mean, as expected for (the square of) a gaussian random variable. The intermediate case is discussed in Scott et 
al. 1993. 
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and so fall off rapidly with I. Each y\ is an independent random variable with distribution 

1-1/2 -y 

with 6 the Heaviside (or step) function. The correlation function is thus a random variable with 
distribution given by the convolution 



P(C) = lim ( J] cT 1 ) P 2 o P 3 o • • • o P N (C) 

N ^°° Vr=2 / 



where Pi(x) = P^(x/ci). To simplify this we take the Fourier transform (compare with Cayon et 
al. 1991) 

iv ^°°Z=2 C ' 1=2 

This can be easily inverted. Writing log {^F[P}) = —a^ + ib^ we need to calculate 

P(C) = r duo e- a " cos % - uoC) 

2,71 J -oo 

where 

OO 

a w = ^(/ + l/2) logyl + ^Q) 2 

Z=2 

and 

oo 

by = y^(Z + 1/2) arctan(cJQ) 

Z=2 

This expression is exact and the integral can be efficiently done numerically. One only needs to find 
a w and once for each set of q to evaluate the distribution for all C. The results shown in Figure 
1 come from taking 9 = and (af) = (/ + 1/2) -1 . To include additional (independent) components 
of the spectrum one just has a w = J2i a$ an d ^ = J2i b$ where a®, are computed from cf^ as 
above. 
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/ 


% Uncertainty 


10° 


7° 


3.8° 


2 


63 


1.00 


1.00 


1.00 


3 


53 


0.69 


0.70 


0.71 


4 


47 


0.51 


0.53 


0.55 


5 


43 


0.40 


0.43 


0.45 


6 


39 


0.32 


0.35 


0.37 


7 


37 


0.25 


0.29 


0.32 


8 


34 


0.20 


0.25 


0.28 


9 


32 


0.17 


0.21 


0.25 


10 


31 


0.13 


0.18 


0.22 


11 


29 


0.11 


0.15 


0.20 


12 


28 


0.09 


0.13 


0.18 


13 


27 


0.07 


0.12 


0.16 


14 


26 


0.06 


0.10 


0.15 


15 


25 


0.04 


0.09 


0.13 


2^1=2 




4.19 


5.01 


6.48 



Table 1: The relative uncertainty and expected size of the multipoles measured by the correlation 
function from COBE (smoothed on 10° and 7° scales) and from MIT (smoothed on 3.8° scales) for 
an n = 1 Sachs- Wolfe spectrum. For n < 1 the multipoles drop off faster with I. The quadrupoles 
are normalized to 1 and the sum I — 2, ■ ■ ■ , oo is listed below each column. The uncertainty shown 
is the second moment of the x^+i-distribution for af. 
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Figure Captions 



Figure 1 Results of a Monte-Carlo calculation of AT/T for the (a) COBE and (b) MIT experiments 
assuming (af)s — (I + 1/2) -1 and (o; 2 )t = 0. Also shown are the results obtained using the 
method of appendix B and fits to a \ 2 form (dotted line) and a gaussian (dashed line) with 
the same mean and variance as the Monte-Carlo distribution. 

Figure 2 Power spectra k 3 W^(k) for a CDM model with Q = l,h = 0.5, n = 1 and Q B = 0.03 (solid 
line) and 0.10 (dotted line) as a function of I = 2k/ Hq. Also shown are the window functions 
Wi (as described in appendix A) versus multipole number / for COBE, MIT, SP91 and MAX. 

Figure 3 Predictions for SP91, assuming the COBE central value of AT/T, as a function of n. Error 
bars are the COBE cosmic variance induced 90% confidence level uncertainties, neglecting 
experimental errors and the cosmic variance associated with SP91. Points are given for both 
scalar and tensor contributions (solid) and scalar only (dashed). Each group of points is at 
the same n, the offset is for ease of viewing only. Higher points represent Vl B = 0.10 and low 
points VL B = 0.03. 

Figure 4 Predictions for MAX, assuming the COBE central value of AT/T, as a function of n. Error 
bars are the COBE cosmic variance induced 90% confidence level uncertainties, neglecting 
experimental errors and the cosmic variance associated with MAX. Points are given for both 
scalar and tensor contributions (solid) and scalar only (dashed). Each group of points is at 
the same n, the offset is for ease of viewing only. Higher points represent Q B = 0.10 and low 
points £l B = 0.03. 

Figure 5 The probability distribution of the ratio of tensor to scalar contributions to the CMB quadrupole 
anisotropy, compared to the ratio of expectation values for these quantities. Also shown is 
the integrated probability distribution, useful for computing confidence limits. 
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